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Abstract. Collisions between nuclei at ultrarelativistic energies produce a color-deconfined 
plasma that expands explosively and rapidly reverts to the color-confined (hadronic) state. In 
non-central collisions, the zone of hot matter is transversely anisotropic and may be "tilted" 
relative to the direction of the incoming beams. As the matter cools and expands into the 
vacuum, the evolution of the system shape depends sensitively on the dynamical response of 
the plasma under extreme conditions. Two-pion intensity interferometry performed relative to 
the impact parameter can be used to measure the approximate final shape of the system, when 
pions decouple from the system. We use several transport models to illustrate the dependence 
of the final shape on the QCD equation of state and late-stage hadronic rescattering. The 
dependence of the final shape on collision energy may reveal non-trivial structures in the QCD 
phase diagram. Indeed, the few measurements published to date show a tantalizing behaviour 
in an energy region under intense experimental and theoretical scrutiny, as signatures of a first- 
order phase transition may appear there. We discuss strong parallels between shape studies in 
heavy ion collisions and those in two other strongly-coupled systems. 
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1. INTRODUCTION 

Color confinement is the most unique and important feature of Quantum Chromodynamics 
(QCD), befieved to be the correct field theory of the Strong Interaction. While the symmetries 
and much of the dynamics of the interaction are understood, a complete understanding of the 
confinement mechanism is hindered by the difficulty of theoretical calculations in a non- 
perturbative regime. Ultra-relativistic heavy ion collisions serve as an ideal, if fleeting, 
laboratory for the study of deconfinement. The system created in such collisions is large 
relative to the hadronic scale; it is characterized by energy densities and temperatures 
sufficient to generate a deconfined multiparticle state. 

Instead of concentrating on the microscopic mechanism of confinement, a fruitful 
strategy has been to focus on "soft" (low momentum or large length scale) observables to 
probe the properties of the bulk matter itself. Of particular interest is the QCD equation of 
state (EoS), quantifying the relationships between intrinsic quantities such as energy density 
and temperature. The nature of these relationships depends on the phase, and as in condensed 
matter physics, much may be learned by focusing on the transitions between phases. 

In nuclear collisions with center-of-mass energy ^/sNN ~ 100 GeV, there is considerable 
evidence that the matter passes into the deconfined phase during some part of its evolution (T] 
El El SI 13. Data do not show evidence of a latent heat associated with the transition (e.g. 
prolonged lifetime ||6| |71 or zero-pressure mixed phase), and there is general agreement 
that the transition is a smooth cross-over at T ^ 150 MeV and jub 30 MeV. Generic 
considerations (H lead to the expectation that a critical point and first-order phase transition 
will be found elsewhere- at lower temperatures and larger chemical potentials- on the T —jub 
plane. Figure [T] shows a schematic of the QCD phase diagram | 9 |. 

Locating and studying these landmarks on the landscape of the QCD phase diagram 
would reveal much about the interaction, including fundamental quantities such as latent 
heats, critical exponents and universality class. Creating a system that samples these more 
interesting conditions requires lowering the collision energy. For this reason, experiments at 
RHIC have embarked on a major program to map out the energy dependence of experimental 
observables, under fixed detector and analysis conditions |9|. Ideally, this exercise parallels 
that of a condensed matter physicist, precisely controlling the temperature of a material and 
measuring its resistance, a precipitous drop clearly marking the transition to a superconductor. 

However, unlike the condensed matter lab sample, the system created in RHIC collisions 
is highly dynamic and far from infinite; phase transition signals in heavy ion collisions will be 
more subtle. The dynamics of any bulk system are dictated by its EoS, which encodes the non- 
trivial landmarks on the phase diagram discussed above. Thus, it is particularly interesting to 
identify bulk observables sensitive to the dynamics and EoS. 

Two-particle intensity interferometry probes the final- state geometry of a nuclear 
collision at the femtometer scale 1 10 |. Detailed measurements of this geometry as a function 
of transverse momentum (pr), rapidity (y) and particle mass, have revealed a rich spacial 
substructure of the system due to pressure-driven bulk collective flow. It has long been 
recognized that azimuthally- sensitive studies in non-central collisions generally yield more 
insight than do azimuthally-integrated measurements. Unfortunately, while momentum- space 
measurements of azimuthal anistropies are plentiful at several collision energies, there are few 
such measurements in coordinate- space. 

In this paper, we compare anisotropic shape measurements in heavy ion collisions with 
similar studies in two other strongly-interacting systems at very different spatial and thermal 
scales. We discuss the physics associated with anisotropic shapes in these collisions, and 
how these shapes are measured. Several theoretical transport calculations are used to show 
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Quark-Gluon Plasma 




Figure 1. Sketch of the QCD phase diagram in the T — /ub plane. An estimate of the critical 
point and first-order transition line is indicated, as are possible trajectories of systems created 
in collisions with energy y/sf^f^. From |9jj. 



the sensitivity of the final shape of the source on the underlying physics, and are compared 
with existing measurements. Particularly interesting is an "anomalous" shape measurement 
at about s/s^n ^ 20 GeV, an energy region where several threshold-like behaviours have been 
reported ifTTIl . Given the extreme paucity of shape measurements, however, it is difficult to 
conclude much from the data at this point. Rather, our work represents a call for a much more 
detailed shape analysis at low RHIC energies. 

2. Anisotropic shapes as a probe of three strongly-coupled systems 

As evidenced by this special volume, there is increasing recognition of connections between 
the physics of ultra-relativistic heavy ion collisions and other quite different fields. In this 
section, we briefly discuss three systems of vastly different scales and decidedly distinct 
physical constituents; these are listed in table[T] Despite their differences, they share a striking 
resemblence. In all cases, a strongly-interacting system is initially prepared in a spatially 
anisotropic state and then allowed to evolve. The shape at a later time reveals important 
physics driving the dynamics of the matter. 

2.7. Anisotropic shape evolution in a cold atomic gas 

The first connection to the bulk evolution in heavy ion collisions and that of cold atomic 
systems, was pointed out several years ago 1 12^ J3J. In particular, in measurements by O'Hara 
et al a degenerate Fermi gas of ultra-cold ^Li atoms is held in a spatially-anisotropic 
magnetic trap. The trap is removed, and the shape of the system is measured at a later time; 
Figure [2] shows the state of the system at different times after the trap is released. (The 
measurement of the system shape actually destroys the cold gas, so the panels of fig. |2] are in 
reality different gas samples released from identical traps.) 

For the system depicted in Figure |2j a pumping laser has been used to maximize the 
inter-atom interaction cross section. Thus, like the ultra-hot partonic system created at RHIC, 
the ultra-cold degenerate gas is a strongly-coupled quantum fluid, starting from an elongated 
ellipsoidal configuration and allowed to expand into the surrounding vacuum. As with the 
QGP at RHIC, the higher pressure gradients along the short direction of the initial shape 
(horizontal in the figure) lead to a stronger expansion in that direction; with the passage of 
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system 


T(K) 


length (m) 


time (s) 


cold atoms 


10-^ 


10-4 


10-^ 


electrical plasmas in crystals 




10-7 


10-12 


heavy ion collisions 


1012 


10-15 


10-24 



Table 1. Characteristic scales for three strongly-coupled systems. Despite their widely 
different scales and matter characteristics, their evolution and methods used to study them 
are strikingly similar. They are prepared in an anisotropic state, expand hydrodynamically, 
and their final shape is studied to determine their underlying physical properties. 

time, the system becomes more round and after some time (~ 600 ^s in the figure) even 
reverses the sense of its elongation. 

The similarity between these measurements and those of "elliptic flow" in heavy ion 
collisions [15| has been noted by others |12, 13 1. However, measurements of elliptic flow 
are restricted to anisotropics in momentum space; a more direct connection is made in space- 
time. Eleven orders of magnitude smaller and 21 orders faster, the heavy-ion analog is shown 
in the hydrodynamical calculations of Figure|2]| 12 |. Here, the initial anisotropy of the system 
is generated by the finite impact parameter of the collision. As with the cold atomic gas, the 
system rapidly expands preferentially along its shorter axis. 

In both systems, the strength of the expansion is driven by pressure gradients which are, 
in turn, determined by the energy density through the equation of state and thermodynamic 
state of the system. The inversion of the aspect ratio in the cold gas system, seen about 
700 JUS after its release, signals a strongly interacting phase, semi-quantitatively understood 
as a superfluid state [14 , 16] . For the heavy ion case, the calculations on the left and right in 
Figure |2] begin with identical initial energy distributions; only the EoS is different. If the EoS 
of a massless gas is assumed (right column), the pressure is large and the expansion rapid. For 
an EoS featuring a first-order phase transition, the pressure gradients are initially large (in the 
QGP phase), then very small (as the system passes through the mixed phase) and finally of 
moderate strength (in the confined phase). Since different regions of the system pass through 
these phases at different times, the flow pattern is complex. What is clear is that the freeze- 
out shape in coordinate space is very sensitive both to the EoS of the hot matter of interest, 
and to the timescale over which its evolution takes place. It is this shape that is extracted by 
azimuthally-dependent femtoscopic measurements discussed in this paper. 

2.2. Anisotropic shape evolution in an electrically-deconfined plasma 

Recently, there has been much activity in the study of dynamic condensed matter systems 
under extreme conditions. This work bears even greater similarity to our study of spatial 
anisotropics in heavy ion collisions. 

The study of condensed matter systems under extreme pressure has a long history; 
recently, using diamond anvils, pressures of up to 0.1 TPa (10^ atm) can be achieved under 
static conditions in the laboratory, resulting in a surprising diversity of new materials ifTSl . 
Exploring even more extreme conditions- 10 TPa or larger- requires explosive generation 
of a transient system such as are done at the National Ignition Facility 1 19] or in table-top 
experiments in which sub-ps laser pulses generate "microexplosions" under the surface of 
Sapphire crystals 1171 crystals or fused silica ll20l . 

The study of these microexplosions parallels strikingly the study of femtoexplosions 
in heavy ion collisions. In the initial state, the matter is in the charge-confined (atomic) 
state. Upon rapid deposition of extreme energy density (10^^ J/m^), a charge-deconfined 
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Figure 2. Time evolution of spatial anisotropy two strongly-coupled systems. Left: A 
degenerate Fermi gas of ultra-cold Li atoms released from an anisotropic trap. From 1 14l. 
Reprinted with permission from AAAS. Right: Hydrodynamical calculation of the evolution 
of a Au-i-Au collision at ^/sNN = 130 GeV. Evolution on the right corresponds to an equation 
of state (EoS) for an ideal massless gas. On the left, the EoS includes a first-order transition 
between hadronic and QGP phases. From 1 12 1, reprinted with permission. 



plasma is generated within a few fs, at temperatures of 10^ — 10^ K. The plasma expands 
rapidly (~ ps), cooling as it does so, and returns to charge-confined degrees of freedom. 
Plasma hydrodynamics and two-component "blast- wave" pictures |21 1 are used to describe 
and understand the source evolution 1171 . 

With huge changes in physical scales and "color charge" replacing "electric charge," 
the above describes the situation with RHIC collisions rather well, down to the blast-wave 
parameterizations |22|. In both cases, too, the final-state anisotropy carries important physical 
information. The anisotropic final-state geometry of a microexplosion is measured directly 
by a scanning electron microscope; c.f. Figure [3] In a heavy ion experiment, it is the final- 
state momenta that are directly measured, and azimuthally- sensitive two-particle intensity 
interferometry must be used to measure the coordinate- space geometry. 

Since the first proof-of-principle microexplosion experiments, there has been 
considerable activity to extract the equation of state of the matter- the plasma state, phase 
transitions, etc. The approach taken is essentially identical to the one we now propose at 
RHIC: to measure the final-state anisotropy as the initial energy of the system is varied, and 
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Figure 3. Anisotropic final state of the evolution of dynamic plasmas. Left: Scanning electron 
microscope images of the aftermath of a microexplosion in sapphire produced by a 150-fs, 
100-nJ laser pulse. Reprinted with permission from 1 17 1; copyright (2006) by the American 
Physical Society. Right: Simplified parameterization of the freezeout distribution in a heavy 
ion collision, as an ellipsoid tilted with respect to the beam axis. Inset: Projection of the 
distribution in the transverse plane. 



compare the results to transport calculations with different EoS. Even the cryptic names of the 
equations of state 1231 (QEOS, SESAME 7387, etc) are reminiscent of those used in RHIC 
studies. 



2.3. Anisotropic shape evolution in hot QCD matter 

The case on which we shall focus henceforth is the anisotropic evolution of the hot matter 
generated in the overlap zone of two colliding heavy nuclei; this is indicated in the right 
panels of Figure |2] Here, we introduce the anisotropics of interest and the physics driving 
their evolution. The situation with heavy ion collisions bears more resemblence to that of the 



microexplosions of section 2.2 than to the cold atoms discussed in 2.1 since the experimenter 
cannot freely choose the time to measure the system anisotropy. When particles decouple 
from the medium created in a heavy ion collision, they are said to "freeze out." Only the final 
state of the system- after it has expanded and frozen out- is available for examination; its 
temporal evolution must be modeled. 

The anisotropy of the hot zone in a heavy ion collision has two sources. Firstly, the 
beam direction (z) is clearly special; both in momentum- and coordinate- space, the hot source 
is extended in the z. Collisions at finite impact parameter break the remaining symmetry in 
the azimuthal variable around the beam direction. The so-called reaction plane is the plane 
spanned by the impact parameter (oriented in the x-direction in this work) and the beam 
direction. Figure |3] shows a plausible if simplistic sketch of the hot matter produced in a 
non-central heavy ion collision, containing the minimal set of possible anisotropics- different 
length scales in each direction, and a tilt of the source away from the beam axis. 

O f par ticular interest is the transverse eccentricity of the source, mentioned already in 
section |2.l| This eccentricity may be quantified by 8 = (a^ — c>x) / + where a^^ 
are characteristic scales of the system in and out of the reaction plane, respectively, and will 
be discussed in more detail shortly. As discussed there, and seen in Figure |2j the final state 
eccentricity is determined by both the anisotropic pressure gradient and the system lifetime; 
increasing either or both of these results in a lower (possibly negative) 8. 

The other major feature of the freezeout distribution is the tilt of its major axis, relative to 
the beam direction. Such tilts are ubiquitously produced in three-dimensional simulations of 
heavy ion collisions. At low energies (-^/snn ^ 4 GeV), 0^ ^ 30° |24|; its sign discriminated 
between competing explanations of momentum- space anisotropics for charged pions L25 J . 
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Figure 4. Important aspects of a tilted QGP. Left: The tilt of the bulk "twisted sQGP" in 
general will not coincide with that of the hard partonic collisions. The interplay between the 
two can distinguish different initial-state scenarios and strongly affects jet quenching signals. 
Reprinted with permission from |26|; copyright (2006) by the American Physical Society. 
Right: The "anti-" (or third) flow component (blue arrows) arises from the expansion of the 
tilted source along its short axis and partially cancels the normal flow (red arrows). Reprinted 
with permission from l30l : copyright (2000) by the American Physical Society. 



No tilt meaurements have yet been made at ultra-relativistic energies, but several theorists 
have pointed to its importance. In what they have termed the "twisted sQGP," Adil and 
collaborators (261 123 emphasize the importance of the tilt of the strongly-coupled quark- 
gluon plasma created in non-central collisions; c.f. Figure]?] In particular, the source of hard 
partonic scatterings (leading to jets) will have a different tilt than that of the plasma itself. The 
interplay between the two tilts is crucial to obtain fully three-dimensional jet tomography f2Jl. 
The hard-partonic tilt depends on the initial-state model, while the tilt (or "twist") of the sQGP 
can be measured by techniques discussed in the next section. 

The tilt of the sQGP leads to important signals in the bulk sector, as well. In three- 
dimensional fluid dynamic calculations, a tilted hot zone generates a collective structure 
known as the third flow |i28) |29| or anti-flow |30|, as it expands preferentially along its 
shortened axis; c.f. Figure]?] This tilt or "torque" can arise naturally in a wounded nucleon 
initial condition 1311 1321 . This third flow component combines nontrivially with "normal" 
directed flow, leading to partial cancelation at low energies 1301 as suggested in the right 
panel of Figure]?] while the third flow dominates at the LHC f33l. Directed, or first-order flow 
signals in momentum space ifTSl are among the most important bulk phenomena, sensitive to 
the earliest, densest stages of the collision 1 34] which may or may not be thermalized. Spatial 
tilt measurements- even if only of the freezeout distribution- will need to be combined with 
momentum-only analyses to disentangle the dynamics of this crucial stage of the collision. 

Generic expectations for the collision energy dependence of freezeout shapes seem 
straightforward. The eccentricity, £, is affected by pressure and timescale. One expects both 
the lifetime and the energy density of the system to increase with increasing ^/s^n- Thus- if 
the relationship between pressure and energy density (the EoS) remains fixed- it is natural to 
expect 8 to decrease monotonically with ^/snn- The tilt is a manifestly non-boost-invariant 
aspect of the QGP created in the collision. Directed flow measurements at all energies confirm 
that the dynamics of heavy ion collisions are never, strictly speaking, boost invariant. Even the 
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hope that the system is "essentially" boost-invariant at midrapidity may be easily shattered if 
a finite tilt angle is measured there. Nevertheless, due to the increased elongation of dynamics 
along the beam direction, it is natural to expect a monotonic decrease of 0^ with -y^s^N as well. 



3. Measuring source anisotropy in heavy ion collisions 

In the two parallel cases discussed in section [2j the final spatial source anisotropy is measured 
directly. Of course, the spatial and temporal scales involved in a heavy ion collision render 
such measurements impossible. Instead, spatial sizes and shapes are extracted via two-particle 
femtoscopy 1 10 |. This technique exploits the connection between the measured two-particle 
relative momentum correlation function C{q) and the spatial separation distribution S{f), 
according to the Koonin-Pratt equation [ 35 1, 



C{q) = ld'r'S{r'mq\r')f 



(1) 



where (|)(^,r) is the two-particle wavefunction as a function of the relative momentum and 
separation in the pair center of mass system, q = pi—p2f'= x\ —x'2. Femtoscopy has been 
used extensively to map the space-time structure of heavy ion collisions for a quarter century; 
details of the method and the physics learned have been discussed elsewhere, e.g. 1 10|. 

Significantly complicating any femtoscopic analysis is the fact that, in Equation [T} C{q) 
and Sir) may- and in reality always do- depend on the pair momentum K=^{p\^p2). This 
is a consequence of space-momentum correlations that arise from collective flow and means 
that particles emitted with a given velocity measure only part of the source. Reconstructing 
the "whole" source from the various fragments is highly non-trivial and ultimately model- 
dependent. For the moment, we ignore this fact and assume that the pion emission probability 
from any space-time point of the source is independent of its momentum. We return to the 
issue of position-momentum correlations at the end of the section. 

The most common femtoscopic analysis in heavy ion collisions correlates identical 
pions. Often, the separation distribution is assumed Gaussian; in this case, ignoring Coulomb 
and other complications 1 10], the correlation function itself is Gaussian, and fitted with 

C(^) = l+^exp[- ^ qiqjRlX (2) 

The indices indicate the components of the relative momentum vector in the Bertsch-Pratt 
"out-side-long" coordinate system, in which the "long" direction is parallel to the colliding 
beams, "out" is parallel to the pair transverse momentum Kt = \{pt,i + Pt.i), and "side" 
is perpendicular to "out" and "long." Use of this coordinate system is motivated by the fact 
that spatial and temporal aspects of the source are more easily disentangled |61 1361 . The 
fit parameters Rfj are squared "HBT radii" that characterize the three-dimensional size and 
shape of the separation distribution. 

In an azimuthally sensitive femtoscopic study, these HBT radii are measured as a 

function of the pair angle = Z (^Kj^b^. Since the "out" and "side" directions rotated 
relative to the x and y directions by even in the simplest case of Figure [3j one 
expects oscillations in the Rfj{(^p). Indeed, such oscillations in non-central collisions 
have been clearly observed; Figure |5] shows HBT radii measured at the lowest (^/sNN = 
2.2 GeV) |24| and highest (-s/snn = 200 GeV) |37| energy collisions explored by these 
analyses. Suppression of oscillations due to the finite reaction-plane resolution may be 



Shape analysis of strongly-interacting systems: the heavy ion case 



9 



Au+Au E(beam)=2 AGeV 




Figure 5. Oscillating HBT radii at the lowest and highest measured energies. Left: ^/sNN = 
2.35 GeV Au+Au collisions with impact parameter b ^ 5 fm. 1 24 1 Right: ^/s^n = 200 GeV 
Au+Au collisions of varying centrality. Reprinted with permission from 1 37 1; copyright (2004) 
by the American Physical Society. 



corrected for fSE\. Similar to the analysis of momentum flow, Fourier moments of the 
oscillating radii are extracted. 

^l,n = { ( ^2(^p)cos(^(^p))(a/ = o,s,l,ol) {Rl{(^p)sm{n^p)){jLi = os,sl).(3) 

Very few femtoscopic studies to date have been performed relative to the event-wise 
reaction plane; instead, the direction of the impact parameter is ignored, and the resulting 
length scales represent an azimuthal average of all collisions. In this case, the pair emission 
angle (^p is meaningless, and no oscillations are measured. Furthermore, the "cross-term" 
radii R^j^i vanish by symmetry f38 1 in these analyses. 

A Gaussian functional form is the most commonly-used parameterization of the pion- 
emitting source in heavy ion collisions; it is characterized by three scales: a lifetime, and 
spatial scales in the longitudinal and transverse directions. The minimal generalization in the 
azimuthally-sensitive case is a Gaussian ellipsoid with three principle axis lengths o^j^^, a 
timescale a^, and a tilt angle 0^ relative to the beam direction. 

/ (xcose^-zsine^)^ / (xsine^+zcose^)^ \ 
/(W,r)~exp(^ ^ (4) 

The primes on ay and c^f denote that these are the lengths of the primary axes of the ellipse, 
the source widths in the tilted coordinate system (c.f. Figure [Sj. 

In the simplest case in which Equation |4] describes the pion-emitting source, its 
anisotropy parameters are simply related to the measured HBT radii. The spatial eccentricity 
along the beam axis is given by 1 22 1 

^y-^l ^?2 

^=-0 ^=2-^. (5) 

+ Rio 

Here, (without the prime) denotes the width of the source in the x-direction, as seen along 
the beam direction, not along the (tilted) major axis of the ellipse. 
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If the first-order oscillations in R^^ is measured, the tilt angle is estimated by combining 
several Fourier coefficients 12511391 , 

and the transverse eccentricity in the "natural" frame tilted relative to the beam axis is 1391 
, - al K2 (1 + ^o^' ^s) + (rIo - Rio) sin' e, - 27?? sin2e, 
^ ^lo (1 + cos2 e,) + [2Rl^ ^Rj^j sin^ 0, + 2/?^ ^ sin2e. 

As suggested in the inset to Figure |3j G^' < Cx for a simple tilted elongated ellipsoid. 
Hence, if the final-state emitting source retains its initial out-of-plane extension (Cy > ay), 
as indicated by measurements so far, the eccentricity measured about the beam axis will be 
smaller than that measured about the tilted axis: 8 < e'. 

At low energies, the direction of the impact parameter b can be estimated easily, thanks to 
a relatively strong first-order anisotropy in momentum-space- the so called "directed flow;" 
in this case, (^p has a meaningful range [0,27l]. At RHIC energies, on the other hand, the 
first-order momentum- space anisotropy is weak, while the second-order momentum- space 
anisotropy ("elliptic flow") is much easier to measure. Thus, at the higher energies, only 
the plane that contains b is defined, but not the direction of b itself; this corresponds, in 
Figure [3j to identifying the yellow reaction plane, but not distinguishing ±x. In this case, (^p 
is measured only modulo K and first order oscillations like R^^ ^ cannot be measured. Spatial 
information on the source tilt 0^ and eccentricity in the source's natural coordinate system 
are inaccessible in such analyses. Similarly, two-dimensional transport calculations such as 
2+1 -dimensional boost-invariant hydrodynamics 1 12] are implicitly blind to any tilt structure. 

As mentioned at the beginning of this section, strong position-momentum correlations, 
due to collective flow or other sources, imply that pion pairs measured at a given do 
not sample the entire pion-emitting zone, but only a selected "homogeneity region" r40l. 
In principle, the correspondence between the homogeneity regions and the "whole" source 
can be almost arbitrary, so that extracting the shape of the latter through measurement of 
the former is neccessarily model dependent. However, studies with reasonable blast- wave 
parameterizations |22 | and realistic transport calculations |39 |- both of which feature strong 
flow and non-trivial correspondence between homogeneity regions and the entire source- 
indicate that Equations [5][7] are good to a model-dependent systematic uncertainty of ~ 30%. 

4. Compilation of experimental results 

Three experiments, listed in Table [2j have published azimuthally- sensitive pion HBT radii. 
All estimated the impact parameter of the collision based on charged particle multiplicity. 



Experiment 


^/s^ (GeV) 


centrality (%) 


rapidity 


AGS/E895 |24J 


2.35, 3.04, 3.61 


(7.4-29.7) 


\y\ < 0.6 


SPS/CERES (m 


17.3 


(7.5-10)0(10-15)0(15-25) 
and (10-15)0(15-25) 


-!<)'< -0.5 


RHIC/STAR 021 


200 


(5-10)0(10-20)0(20-30) 
and (10-20)0(20-30) 


\y\ < 0.5 



Table 2. Measurements of the anisotropic shapes from heavy ion collisions. The third column 
indicates which centrality bins were averaged, to obtain the shape parameters of Figures [6] 
andlT] See text for details. 
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Figure 6. Source tilt relative to the beam direction vs. energy for midcentral heavy 
ion collisions. See sections |4] and [5] for a discussion of the experimental data and model 
calculations, respectively. 



Since system anisotropy clearly depends on the impact parameter, it is important to compare 
collisions with similar centrality. In order to best compare results from the 7.4 — 29.7% 
centrality cut of E895 (corresponding to /? = 4 — 8 fm), several centrality cuts were combined, 
for the higher energy measurements. 

It is worthwhile to describe in detail how centrality bins were merged, since the 
comparison between shapes from different collision energies is important for our message. 
The most relevant centrality bins reported by the STAR Collaboration are for 5 — 10%, 
10 — 20% and 20 — 30% of total cross section. We combine data from these three bins as 

1 V pSTAR I 9 V cSTAR i o v pSTAR 
STAR ^ ^ ^ ^(5-10) + ^ ^(10-20) + ^ ^(20-30) 

^(5-10)0(10-20)0(20-30)— l_^2 + 2 * 

Here, the weighting factors (1,2,2) account for the fact that the 5 — 10% bin contains half 
the number of events as either of the other two bins listed. This selection includes more 
central events- in particular, events in the 5 — 7.4% centrality range- than are included in the 
E895 cuts. Therefore, the value £^5^k))0(io-2O)0(2O-3O) ~ 0.081 ±0.006 should be considered 
a lower bound on the shape compared with E895. An upper bound may be obtained by 
combining only the two more peripheral bins: £^1^0^0)0(20-30) ~ 0.094 ± 0.007. 

The relevant centrality ranges reported by CERES are 7.5 — 10%, 10 — 15% and 
15 — 25% of the total cross section. A fair range to use in our comparison is between 
^M)e(io-i5)e(i5-25) = 0.035 ±0.018 and e'^'lTim^s-is) = 0-043 ± 0.020. 

All measurements focused on low-momentum pions (pr ^ 0.25 GeV/c), for which the 
formulae [5][7] work best ||22| . STAR and E895 measurements center on midrapidity, where 
participant contributions should be maximal, while the CERES measurement is somewhat 
backwards in the center of mass frame. Since HBT measurements typically vary slowly with 
rapidity, this difference is unlikely to affect CERES' shape estimation, but a measurement at 
midrapidity would provide a better comparison with the other experiments. 

E895 measured HBT radii relative to the first-order event plane (i.e., the direction of the 
impact parameter); results for s/snn = 2.35 GeV are shown on the left panel of Figure [5] The 
spatial tilt is shown in Figure |6] The tilt is strikingly large at these low energies and drops 
with energy, consistent with the expectation |42 | that collisions become increasingly boost- 
invariant (at least near midrapidity) with increasing energy. It will be important to extend tilt 
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Figure 7. Source eccentricity vs. energy for midcentral heavy ion collisions. See sections [4] 
and[5]for a discussion of the experimental data and model calculations, respectively. 



measurements to higher energies, since a finite 0^ is manifestly "boost- variant," even aty = 0. 
If 0^ is more than a few degrees, boost-invariant models may not be valid and would at least 
require double-checking with true three-dimensional calculations. 

Figure [7] shows the measurements from the experiments listed in Table [2j Filled symbols 
indicate 8, the eccentricity relative to the beam axis (c.f. Eq.[5|, while open symbols indicate 
the eccentricity in the natural frame of the source (Eq. [7]), measured only by E895. For the 
CERES and STAR datapoints, the average of the upper and lower bounds discussed above 
are plotted, with the difference between the bounds and the statistical errorbars added, to 
be conservative. The non-monotonic behaviour of 8 (y^is^viv) is intriguing. As discussed in 
Section 2.3 rather general considerations lead to the expectation of a monotonic decrease of 
8 with energy. The unexpected dip in Figure [7] occurs in the energy region in which phase 
transition "threshold" effects have been reported |11| and around which some speculate that 
heavy ion collisions sample the non-trivial features sketched in Figure[T} c.f |9|. 

To contribute our own speculation, we note that such non-monotonic behaviour could 
arise from one of two effects, both related to a first-order phase transition. Firstly, an extended 
lifetime due to the transition would allow the system to evolve further towards a round 
shape (c.f. Figure |2]), causing a dip just around the threshold energy. Using this simplistic 
scenario to explain the data, the CERES datapoint at ^/sNN GeV lies near the threshold 
energy. Alternatively, we may fix the lifetime and consider effects of the stiffness of the EoS, 
quantified by the speed of sound in the medium, = where P and e are the pressure 
and energy density, respectively. At low energies, the system is in the hadronic phase, and 
^ i ; as ^/sNN increases, pressure gradients increase in proportion to the energy deposited 



6' 

in the system; here, 8 would fall with -^/snn. Near the threshold energy, the system may spend 
much of its time in the mixed phase, for which = 0; here, the system shape would evolve 
little from its initial, large value. As the energy increased still further, the system spends most 
of its time in the deconfined plasma phase, for which ~ 1/3, and 8 (y^^vviv) again falls 
monotonically. Using this second simplistic scenario to explain the data, the STAR datapoint 
at s/snn = 200 GeV lies just above the threshold energy. 

However, such simplistic considerations only serve to stimulate more sophisticated 
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treatment with theoretical tools. We consider such tools in the next section. 
5. Transport model calculations 

Transport models are commonly used to simulate the dynamics of a heavy ion collision. They 
fall into two broad categories- hydrodynamical calculations and Boltzmann transport models. 

As the name implies, hydrodyanmic calculations model the nuclear matter in terms of 
one or more fluids, each fluid cell characterized by thermodynamically intensive quantities 
such as temperature, chemical potential, and pressure. Particles per se play no role in the 
dynamics; rather, they appear at the end of the simulation, as the fluid cells produce final-state 
particles according to an approximate prescription B3l . A major advantage of hydrodynamic 
calculations is that the equation of state must be explicitly defined. The relationships between 
bulk quantities- pressure, energy density, temperature, baryon number, etc- close the set 
of hydrodynamic equations, which are otherwise "trivial" continuity and conservation laws. 
These relationships characterize and quantify the state of matter under different conditions 
and are the mathematical implementation of phase diagrams such as that in Figure [T] 

While hydro calculations focus explicitly on bulk quantities, microscopic Boltzmann 
transport calculations take the opposite tack: the creation, scatterings and demise of each 
particle is followed. Such calculations offer several advantages over hydrodynamic ones: 
they are manifestly three-dimensional (important in the present context); assumptions of 
thermalization and equilibrium, justified or not, are not required; event-by-event fluctuations- 
potentially large in such finite systems and of increasing interest in the field- are implicitly 
included; bulk and shear viscosity are automatically included in the calculation; finally, 
the final state of the system- freeze-out- occurs "naturally" on a particle-by -particle basis, 
without recourse to arbitrary criteria (e.g. a fixed local temperature) setting the freeze-out 
hypersurface. The chief drawback of Boltzmann transport calculations is their complexity; the 
dynamics is influenced by the decay and reaction details of several thousand hadron species, 
as well as the string phenomenology needed to describe non-hadronic degrees of freedom in 
the ultra-high density limit at which 2^2 scattering simulations break down. 

We use several transport calculations to simulate the evolution of a heavy-ion collision 
from its initial out-of-plane-deformed shape, to the final shape probed by femtoscopy. 

5.7. Two-dimensional boost-invariant hydrodynamics 

Among the earliest and most exciting successes in the RHIC program was the degree to 
which ideal hydrodynamic calculations predicted the magnitude and details (mass- and pr- 
dependence) of bulk collective behavior in non-central heavy ion collisions |[T2l [TSl 1441 . 
Here, we use one of the most successful models- AZHYDRO | 45 1, a (2+l)-dimensional ideal 
hydrodynamic model. This model uses a common simplifying assumption ||46l |42]| that the 
initial conditions and subsequent dynamics are boost-invariant; all interesting physics takes 
place in the transverse plane. For our purposes, this means that the model assumes a non-tilted 
source (0^ = 0), and only 8 is calculable. 

It is by now rather clear that boost-invariant ideal hydrodynamic calculations are invalid 
at energies below about y^^vviv ^ 20 GeV 1 12 1 as (i) the densities achieved are inconsistent with 
a zero-mean-free-path approximation, and (ii) the system is not boost-invariant. Nevertheless, 
we will perform calculations over a wide range of energies, including very low ones, to map 
the excitation function and explore the EoS dependence of the shape parameters. 

In order to extend the model into both the LHC high energy regime as well as low energy 
heavy-ion collisions, the initialization routine must be tuned to the appropriate collision 
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energies. The initial temperature distribution was parameterized through the initial transverse 
entropy profile geometrically using an optical Glauber calculation discussed in 1121 . 

Within the optical Glauber model calculation, the density of wounded nucleons (A^^) 
and binary nucleon-nucleon collisions (A^bin) are estimated in the transverse plane. The total 
entropy density is then a superposition of the "soft" wounded nucleon density and the "hard" 
binary collision density, appropriately scaled f4T, "48 1 to match both the charged hadron 
multiplicity and centrality dependence observed in experiments (25% hard contribution). 
The multiplicity is simply the momentum integrated particle distribution at a given rapidity 
y (typically at midrapidity, y = 0) |12j i47J. Furthermore, the for consistency, as the 
initial entropy density (and consequently the initial temperature) changes for each collision, 
the thermalization time Xq is correspondingly changed in order to keep the "uncertainty 
relationship" (Xq^o ^1) constant fT2^,'49 |. The multiplicity is then matched to the expected 
charged multiplicity of each collision energy. 

To generate the observable momentum distribution of particles, a freeze-out criterion 
is needed at the end of the hydrodynamic simulation. The Cooper-Frye formalism |43 | 
postulates a sudden transition from a perfect local equilibrium to free- streaming for all 
strongly interacting particles in a particular fluid cell at a given kinetic freeze-out condition. 
While this can be described via a dynamic comparison of scattering rates, here we use a 
constant energy density ^ = 0.075 GeV/fm^ which corresponds to a constant temperature 
isotherm which varies depending on the equation of state used to describe the matter. For all 
equations of state used in these simulations, the freeze-out temperature is Tf ^ 130 MeV. 

In order to study a range of possible phases of matter, three distinct equations of state 
were used in the hydrodynamic simulations. These are discussed in |45| and represent a 
purely hadronic state of matter (EOS H), a pure quark/gluon gas (EOS I), and a Maxwellian- 
constructed EOS Q which contains a first order phase transition between EOS I and EOS H. 
In B31 , EOS Q was used exclusively to describe the matter created at s/s^n = 130 GeV. 

HBT radii were extracted by fitting projections of two-pion correlation functions, 
calculated according to the method of |50|. These radii were then used in Equation |5] to 
extract the final-state eccentricities shown in Figure |7] For each EoS used, the eccentricity 
monotonically decreases as a function of energy, an effect both of increased system lifetime 
and pressure as ^/sNN increases, as discussed in section [23| 

There is considerable sensitivity to the EoS used in the calculation. Using the 
initialization procedure discussed above, use of the stiff equation of state, EoS-I with c^ = ^, 
results in a much more out-of-plane shape- i.e. one that has not evolved much from the initial 
overlap shape. The shape evolves considerably more when using the softer EoS-H (cj = ^). 
Since, for a given energy density pressure gradients are proportional to c^, these results suggest 
that effects of system lifetime dominate over those of pressure, in these calculations. This 
conclusion is consistent with the results when using EoS-Q. These shapes track closely with 
those of EoS-H for low ^/sNN where the system is dominated by the hadronic phase. At around 
^/sNN^ the threshold effect of the "soft" mixed phase (cj = 0) become apparent, increasing the 
system lifetime and further decreasing 8. 

It is clear from Figure [7] that the phase transition is needed to explain the single RHIC 
datapoint- at least in this model. Extension of the shape excitation function to higher (LHC) 
energies will be important to further constrain the EoS. However, many of the most important 
features of the EoS may be manifest at lower energies, such as those currently being explored 
in the RHIC energy scan. The validity of (2+l)-dimensional models breaks down at these 
energies, limiting their utility in constraining the EoS with low-energy data. 
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5.2. Microscopic Boltzmann transport 

In this work, we use two related Bolzmann transport codes. The first is the Ultra- 
Relativistic Quantum Molecular Dynamics Model (UrQMD 3.3)211 ES. which includes 
details of physical processes relevant over a huge range of energies. For this reason, it 
is particularly attractive for the present study. UrQMD is a covariant transport approach 
to simulate the interactions between hadrons and nuclei up to relativistic energies. It is 
based on the propagation of nucleons and mesons accompanied by string degrees of freedom 
with interaction probabilities according to measured and calculated cross sections for the 
elementary reactions. Hard scatterings with large momentum transfer are treated via the 
PYTHIA model 1531 . For detailed comparisons of this version to experimental data, the 
reader is referred to ll54ll . 

An earlier incarnation of this model, Relativistic Quantum Molecular Dynamics 
(RQMD) 1^1 was widely used over the more limited energy range of the AGS and SPS- 
^/snn ^3 — 20 GeV. It is satisfying to see good consistency between the older and newer 
versions of the model, in terms of predicted shapes. 

For the UrQMD calculations, 8, and 0^ were extracted by directly fitting the freezeout 
distribution with the functional form of Equation |4] For the older RQMD calculations, the 
model output was processed through Equation[T]to generate correlation functions, which were 
fitted with Equation [2] to extract HBT radii. These radii were then used to calculate shape 
parameters according to Equations [5]|7] For these models, the shape parameters extracted 
using these two methods should be consistent to ^ 30% 1391 , 

The simulations reproduce the very large tilt angles measured at low ^/snn and predict 
a sharp fall-off with energy. The RQMD model features the possibility to include the effects 
of a medium-induced mean field on the trajectories of the hadrons during the collision; 0^ 
is significantly sensitive to this mean-field effect. In particular, the spatial tilt measurements 
are best described when effects of the mean field are ignored ("cascade mode"). This is 
interesting in light of the fact that reproducing the momentum-spsice tilt ("directed flow" or 
vi 1 15 1) demands inclusion of mean-field effects |56|. While measurement of spatial shapes 
already constrain the EoS of hot nuclear matter, combining both coordinate- and momentum- 
space shapes place even stricter constraints on the dynamics. 

The final- state eccentricity, plotted in Figure [Tj reproduces the large 8 (and e') values 
measured at the AGS, with little dependence on the nuclear mean field. At SPS energies 
(s/snn ^17 GeV), tilt angles on the order of 10° are predicted by UrQMD. This is just about 
the point at which the effect of the tilt on the measured eccentricity (c.f. inset of Figure [3]) 
vanishes; i.e. 8^8' for s/snn > 17 GeV in this model. The monotonic decrease predicted 
by the model is not nearly strong enough to reproduce the CERES measurement, but falls 
rather smoothly to closely approach the shape measured at top RHIC energy. At still higher 
energies (e.g. LHC), the model predicts a continued out-of-plane final eccentricity with little 
y^isvviv -dependence . 

5.3. Hybrid hydro dynamic -Boltzmann calculation 

Combined micro scopic+macroscopic approaches are among the most successful ideas for 
the modeling of the bulk properties of HICs (571 ES ESI- The approach that we are using 
here has recently been developed and is based on the UrQMD hadronic transport approach 
including a (3+l)-dimensional one fluid ideal hydrodynamic evolution |[60l|6T| for the hot 
and dense stage of the reaction (621 [63l . To mimic experimental conditions as realistically 
as possible the initial conditions and the final hadronic freeze-out are calculated using the 
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UrQMD approach. The non-equilibrium dynamics in the very early stage of the colHsion 
and the final state interactions are properly taken into account on an event-by-event-basis. 
Furthermore, the hybrid model allows for a dynamical coupling between hydrodynamics 
and transport calculation in such a way that one can compare calculations with various EoS 
during the hydrodynamic evolution and with the pure cascade calculations within the same 
framework. 

The coupling between the UrQMD initial state and the hydrodynamical evolution 
proceeds when the two Lorentz-contracted nuclei have passed through each other, ^tart = 
IRj — 1 1641 . After that, a full (3+1) dimensional ideal hydrodynamic evolution is 
performed using the SHASTA algorithm [60, 61 1. Taking into account all three spatial 
dimensions explicitly in the evolution is important to be able to study the angular dependence 
of HBT radii, since effects like the longitudinal tilt of the event plane can only be consistently 
considered in this way 1391 . 

Serving as an input for the hydrodynamical calculation the EoS strongly influences the 
dynamics of an expanding system. Two different equations of state are used to exemplify the 
differences on the extracted HBT radii due to this external input. One is a hadron gas equation 
of state (HG) with the same degrees of freedom as in the UrQMD approach 1651 . The other 
one is a bag model equation of state (BM) including a strong first order phase transition to 
the quark gluon plasma with a large latent heat 1611 . To see if fluctuations in the initial 
state affect the result differently for different expansion dynamics during the hydrodynamic 
evolution these two extreme cases have been chosen. 

The transition from the hydrodynamic evolution to the transport approach when the 
matter is diluted in the late stage is treated as a gradual transition on an approximated iso- 
eigentime hyper- surface (see (661 1221 for details). The final rescatterings and resonance 
decays are taken into account in the hadronic cascade. 



As with the UrQMD calculations discussed in section 5.2 shape parameters were 



extracted from a direct fit of the freezeout distribution with Equation [4] Large tilts are 
again predicted at low collision energies, with significant sensitivity to the EoS used in the 
calculation. The effect of the first-order phase transition (Hydro[BM]) is clear: while at yjs^ 
even the earliest dense phase of the collision is below threshold to be affected by the phase 
transition. However, at larger energies, the mixed phase reduces the sideward pressure very 
early in the system evolution, reducing 0^. It would be very interesting indeed to measure tilt 
angles at SPS energies. 

That shapes from Hydro [HG]+UrQMD are not identical to those from "pure" UrQMD is 
at first puzzling, given that the EoS in the hydrodynamic phase is that used in the Boltzmann 
model. This is likely a technical issue- particles that would be emitted early in the "pure" 
UrQMD simulation are generally absorbed into the hydro phase in the hybrid model, only to 
reappear at the iso-eigentime hyper- surface mentioned above. The details of this discrepancy 
are still under investigation. Thus, for now it is best not to compare hybrid calculations to 
pure Boltzmann simulations, but to compare one hybrid calculation to the other, concluding 
that the EoS sensitively affects the final-state shape of a heavy ion collision. 

Due to this increased system lifetime, the system eccentricities in the hybrid calculations 
are allowed to evolve to much lower values than those predicted by UrQMD. Indeed, the 
system in its natural rotated frame is essentially round transversely (e' ^ 0), so the large tilt 
can even produce 8 < 0- an in-plane extended source, measured about the beam axis; 8 then 
grows with y^^svviv because 0^ decreases. The eccentricity at SPS energy essentially reproduces 
the CERES result, probably an artifact of the extended lifetime effect discussed above. 
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6. Discussion and Summary 

Locating and studying non-trivial structures (phase transitions, critical points, etc) in the 
phase diagram offers keen insights on the material under consideration. The equation of state 
encodes these structures into dynamical relationships that determine a system's response and 
evolution from an externally imposed initial state. A major program is currently underway 
at the RHIC facility, to vary the collision energy of heavy ion collisions with an eye for non- 
trivial energy dependence of bulk observables that might signal the presence of structures on 
the phase diagram. Several sensitive observables have been proposed and are under study | 9J. 
We have discussed one such observable here- the final-state anisotropic shape of the system 
in coordinate-space. 

The study of such shapes is strikingly similar to similar studies of strongly-coupled 
systems at vastly different scales. The parallel between the evolution of a non-central heavy 
ion collision and that of a strongly-coupled gas of cold atoms released from an anisotropic 
trap, has been noted before, though more in relation to mometum- space, rather than 
coordinate-space, shapes. However, the heavy ion situation bears much more resemblance 
to the study of "micro-explosions" induced by femtosecond laser pulses on crystals. In 
both cases, a charged-confined system is, on very short timescales, raised to an energy 
density sufficiently high to generate a charge-deconfined plasma. The plasma responds 
hydrodynamically- expanding and cooling until the system returns to its original, charged- 
confined phase. The equation of state is extracted by comparing the anisotropic shape of the 
final state with transport calculations with different EoS. The energy scan program at RHIC 
is following the lead of studies of these micro-explosions, varying the initial energy of the 
system as a sensitive way to probe the EoS. 

Whereas the shapes of cold atomic gases or micro-explosions can be measured directly, 
two-particle intensity interferometry is the most direct probe of the space-time structure of 
evolving matter on the femtometer scale. We have discussed how system shapes are obtained 
from such measurements and the very few measurements that have been made to date. 

The shape excitation function, sparse though it is, features an unexpected minimum 
at s/snn =17 GeV, an energy around which other "anomalous" behaviour has been 
reported 1 1 1 1. Based on qualitative arguments, we speculated that a first-order phase transition 
might cause the minimum. This minimum is not predicted by any transport model, even those 
that include a first-order phase transition. It is good to keep in mind, however, that none of the 
models we considered are perfect- the Boltzmann models do not reproduce the large flow seen 
at RHIC, and the two-dimensional hydrodynamic models miss one of the prime anisotropies- 
the source tilt, which is non-vanishing even at SPS energy. In all cases, the sensitivity of the 
final shape of the collision to the equation of state driving the system's evolution, was clear. 

A careful and systematic set of shape measurements at different energies is clearly 
warranted. While some models (e.g. UrQMD and 2D hydro with EoS-Q) predict very similar 
shapes at RHIC (both close to the measured shape), their predictions for higher energies 
diverge strongly. At lower energies, excitation functions of tilts and ellipticities form an 
important part of the energy scan program at RHIC, to search for structures on the fundamental 
phase diagram of QCD. 
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